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Abstract 

I present a simplified model for the gluon Green's function governing high-energy QCD 
dynamics, in arbitrary space-time dimensions. The BFKL integral equation (either with 
or without running coupling) reduces to a second order differential equation that can be 
solved in terms of Bessel and hypergeometric functions. Explicit expressions for the gluon 
density and its anomalous dimension are derived in MS and Qq factorization schemes. This 
analysis illustrates the qualitative features of the QCD gluon density in both factorization 
schemes. In addition, it clarifies the mathematical properties and validates the results 
of the "7-representation" method [1] proposed by M.Ciafaloni and myself for extracting 
resummed next-to-leading-log x anomalous dimensions of phenomenological relevance in 
the two schemes. 
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1 Introduction 



Small-a; resummations in QCD have been extensively investigated in the past years in order to 
improve the fixed order perturbative description of high-energy hard processes in the small-x 
regime, where higher order perturbative corrections grow rapidly due to logarithmically en- 
hanced contributions ~ (og log l/x)". Knowledge of the precise relationship between the fixed 
order approach — based on the coUinear factorization formula and the DGLAP equation [2] 
— and the small-x resummed ones — based on the high-energy factorization formula [3] and 
the BFKL equation [4] — is of course needed for a unified picture of small-x physics, e.g., to 
provide quantitatively accurate predictions in the small-x region, which will be explored by 
next-generation colliders. 

A major aspect of this relationship is the issue of the factorization scheme employed to de- 
fine parton densities and coefficient functions. Fixed order perturbative calculations mostly use 
the (modified) minimal subtraction (MS) scheme in the context of dimensional regularization. 
On the other hand, small-x resummed approaches — being based on fc-factorization [3] which 
involves off-shell intermediate particles with non-vanishing transverse momentum k — are nat- 
urally defined in the so-called Qo-scheme [5], where infra-red (IR) singularities are regularized 
by an off-shell probe whose non-vanishing virtuality Qq plays the role of an IR cutoff. 

The basic relations for the MS ^ Qo scheme change of anomalous dimensions and coefficient 
functions were obtained some time ago [6,7] at relative leading- log x (LLx) order, then improved 
to include next-to-leading-log x (NLx) running coupling corrections [8] and recently extended 
by M.C and myself at full NLx level [1]. The main tool of our analysis [1] was the generalization 
to 4 + 26 dimensions of the 7-representation of the gluon density — a Mellin representation of 
the BFKL solution in which 7 is conjugate to t = log(fc //i^). While for e = the running- 
coupling BFKL equation is a differential equation in 7, for e 7^ it becomes a finite-difference 
equation, whose solution, however, is not unambiguously determined and has been computed 
by using sometimes rather formal manipulations. Despite the sensible physical meaning of the 
procedure and of its results, from a mathematical point of view some steps of our method are 
not fully proven. It is therefore desirable to have at least an explicit example that could confirm 
our method of solution of the finite-difference equation, especially in view of its application to 
compute the anomalous dimensions at full NLx accuracy. 

The purpose of the present work is to devise a non-trivial, physically motivated and solvable 
model which: by providing explicit solutions, illustrates the main qualitative features of the 
real QCD case; 2) can clarify the less understood aspects of the procedure developed in [1] 
and verify the correctness of its results. The model I am going to present is a generalization 
to arbitrary D = A + 2e space-time dimensions of the collinear model [9] used in the past to 
study the interplay between perturbative and non-perturbative QCD dynamics at high energies. 
Starting from the formulation of the LLx BFKL equation in D dimensions, only the coUinearly 
enhanced (fc ^ k' and k -C k') contributions of the integral kernel K{k, k') are kept. Despite 
its poor phenomenological accuracy, this model contains most of the qualitative features of the 
real theory: it is symmetric in the gluon exchange k ^ k' , it generates collinear singularities 
in the e — >• limit, it correctly describes the leading- twist LLx behaviour of the gluon density, 
it includes the running of the coupling. Most importantly, in contrast to the BFKL equation, 
the collinear model can be solved, as a 1-dimensional Schrodinger-like problem. 

Sec. [2] is devoted to the definition of the model in generic number of dimensions. A pre- 
liminary study on the qualitative features of the solution of the master integral equation is 
presented. To this purpose, I briefiy review the two types of running-coupling behaviour that 
are present in 4 -|- 2e dimensions. 

The resolution of the model in the fixed-coupling case is presented in sec. El The ensu- 
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ing integral equation is then recast into a second order differential equation of Bessel type, 
whose solution provides the unintegrated gluon density. The unintegrated gluon density is first 
compared with the known perturbative solution [7] and then used to compute the integrated 
gluon density and anomalous dimension in both the MS-scheme and Qo-scheme. The last part 
of this section concerns the analysis of the Mellin representation of the gluon density and its 
comparison with the corresponding series and integral representations derived in ref. [1]. 

Sec. m includes the one-loop running coupling. In this case the differential equation is solved 
in terms of hypergeometric functions. The analyticity properties of the solution will reveal 
essential in extending the unintegrated gluon density from the IR-free regime — where the 
coupling is bounded — to the ultra-violet (UV)-free regime — where the Landau pole renders 
the integral equation meaningless. The explicit results of the 6-dependent resummed MS and 
Qq anomalous dimensions — which are shown to agree with the known lowest order running 
coupling corrections — provide a strong check for the connection between e-dependence of the 
kernel and 6-dependence of the MS anomalous dimension argued in ref. [1]. 

A final discussion is reported in sec. El 

1.1 Notations 

I distinguish two symbols of asymptotic behaviour: f(x) ~ g{x) means lim^^j-g = k for 
some finite and non-zero k, while f{x) ~ g{x) refers to the special case k = 1. 

The hypergeometric function is denoted by 2-^1 («, b; c; z) = 2-F1 

A citation like [1](2.3) means eq. (2.3) of ref. [1]. 

There are some change of notations between ref. [1] (left side) and this paper (right side): 

^,(fc)e^'^«/r(i + e) m 

-T to. 



2 Formulation of the collinear model 

In this section, I define a simplified model for the gluon density in high-energy QCD with both 
running and frozen coupling constant. After recalling the features of the two running coupling 
regimes, I briefiy discuss the expected qualitative behaviour of the solutions of the model. 



2.1 Motivation of the model 

In high-energy QCD, parton densities and anomalous dimensions are often computed in two 
different factorization schemes, which differ essentially by the regularization of the infra-red 
(IR) singularities. 

• In the so-called Qo scheme [5] , the IR regularization occurs by considering off-shell initial 
partons with non- vanishing virtuality k"^ = —Ql < 0, which plays the role of a momentum 
cut-off; 

• The minimal subtraction (MS) scheme instead, is based on dimensional regularization 
with on-shell initial partons living in Z) = 4 + 25 space-time dimensions, where IR singu- 
larities shows up as poles ~ l/e" and are subtracted from the physical quantities according 
to the MS prescription. 
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The relation between these two schemes can be investigated by including in the defining equa- 
tions for partons both off-shell initial conditions and arbitrary space-time dimensions. 

As for the physical case of 4 space-time dimensions, also in generic D = 4 + 2e dimensions 
the high energy (i.e., small-x) behaviour of cross sections in QCD is governed by the gluon 
Green's function (GGF) Suj,e{k, ko). Here u is the Mellin variable conjugated to x, while k and 
ko are the transverse momenta of the (reggeized) gluons emerging from the impact-factors of the 
external particles [3]. The GGF obeys the integral equation (in the following the dependence 
on the uj variable will always be understood) 

uUk, ko) = 6'+''{k -ko) + J ^^^^ Xeik, k')Seik', ko) (1) 

where the kernel has been determined exactly in the leading-log(x) (LLx) approximation [7] 
and can be conveniently improved to include subleading corrections (in particular the running 
of the coupling). Detailed studies of the ensuing solutions and physical consequences has been 
presented in refs. [7] in the LLx approximation, and in refs. [1, 10] at subleading level. 

It should be noted that the NLa; approximation limits not only the knowledge of the kernel 
but also the method of solution of eq. ([T]). However, it would be desirable to have an 
exact solution of eq. ([T]), even with an approximate kernel, in order to understand the overall 
"non-perturbative" feature of the QCD gluon Green's function. To this purpose, I consider 
a simplified model for %s whose main virtue is to provide a GGF which can be expressed in 
terms of known analytic functions. This toy-kernel resembles the field-theoretical one in the 
coUinear regions k^ -C k''^ and fc^ ^ k'"^, and has already been considered in the past [9] in 
order to study the structure of high-energy QCD dynamics in 4 dimensions. In the following I 
generalize the coUinear model to the dimensional regularized theory, with running coupling as 
well as with fixed coupling constant. 



2.2 Definition of the model 

The coUinear model is defined by the coUinear limit fc^ ^ fc?, of the LLx BFKL high-energy 
evolution kernel 

aC™(fc, k') = ,f ,,,, + Virtual terms ^ = %f\k, k') , (2) 



7r 



(k-ky Tik 



> 



fc< (fc >) being the smallest (biggest) transverse momentum, and g the dimensionful gauge cou- 
pling. By introducing the dimensionless coupling constant a^, the small-x expansion parameter 
a and the logarithmic variable t 

«s = u+e-rn J ^ — ' a = — , t = log — , 3 
(47r)^+^r(l -I- e) TT U /i^ 

we can express both GGF and kernel in terms of the dimensionless quantities fs (the uninte- 
grated gluon density) and defined by 

coSsik, ko) = S'^'%k - ko) + ^^^^m to) (4) 

fc2 Xf\k, k') = ^K'°'\t - t') . (5) 

TT 
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so that one can rewrite eq. ([T]) in the fornix] 

I-+CO 

Mt, to) = ae^* K'°\t - to) + / K'°\t - t')f{t', to) dt' 



ae' 



et 



e(to-t)e*-*o + e(t-to) 



/(t',to)dt'+ / e*-*7(t',to)dt' 



where in the second hne I have substituted the expression of the colhnear kernel 

K'=°"(r) = 0(-r)e" + e(r) , T = t-t' 

stemming from eqs. ([2]) and ([5]). 

A second way to relate this model to QCD is to compare the eigenvalue function 



Xcoii(7) = Xil) 



+00 -1 -1 

e-^"ir^°"(r) dr = - + 

00 



7 1-7 



(6) 
(7) 



(9) 



with the BFKL one Xbfkl = 2'?/'(l) — ip{'y) — ip^l — 7), as in fig. [H Clearly the two eigenvalue 
functions display the same qualitative behaviour in the region around and between the leading- 
twist poles at 7 = 0, 1. 
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Figure 1: Comparison of the colhnear model eigenvalue function (solid-blue) with the BFKL 
one (dashed-red) . 

The colhnear model can be easily generalized to include the running of the coupling. The 
small-x parameter a acquires a t-dependence according to the evolution equation 



= .[a(i) - Bam ■ B - I 



(10) 



where b is the one-loop beta function coefficient (6 = 11/12 — Nf/QNc in QCD). The solution 
of eq. (fTOj) is given by 



at = a(t) 



ae 



et 



1 + a5(e=* - 11 



a(t) 



B 



Ad 



et 



A 



1 - aB 



:ii) 



Note that in dimensional regularization {e 7^ 0) the coupling a(t) has a non-trivial t-dependence 
also in the so-called frozen coupling case corresponding to 6 = 0. Substituting a(t) in place 



^In ref. [1] we adopted a step function Q{t + T), T = —to, instead of K'^°^^{t — to) as inhomogeneous term. 
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Figure 2: Behaviour of the running couphng a{t) in the regular regime aB < 1 (sohd-red) and 
in the Landau regime aB > 1 (dashed-blue) . The straight hne (dotted-green) corresponds to 
the boundary value aB = 1. The case B = is represented by the dash-dotted black curve. 



of ae^* in eq. (Ej), we obtain, after rearranging some terms, the generalization of the colhnear 
model with running coupling: 

fe,b(t, to) 



et 



K^"^Ht, B) = K^"^Hr 



/+00 
K'-\t-t',B)f,4t',to)dt' 
-oo 

BS{t) . 



(12) 
(13) 



2.3 Running coupling regimes 

It is important at this point to realize that the running coupling behaves in two qualitatively 
different ways, according to whether the parameter aB is greater or less than 1. 

• When aB < 1, i.e., 61^ < e/h, the running coupling a{t) is bounded, positive and increases 
monotonically from the IR-stable fixed point a(— 00) = to the UV-stable fixed point 
a(+oo) = 1/-B, as shown in fig. [21 

• When aB > 1, i.e., as > ^/b, the running coupling starts from the positive UV-stable 
fixed point a(+oo) = l/B, then increases and diverges at the Landau point 



--log(-A5) = -log(l-4 

e e \ ab 



(14) 



becomes negative for t < and finally vanishes at t = — 00. This is the situation realizing 
the physical limit e — ^ at fixed b. 

In the former case, the extra-dimension parameter e not only regularizes the IR singularities, 
but avoids also the occurrence of the Landau pole, thus allowing a formulation of the integral 
equation free of singularities. In practice, the strategy of dimensional regularization consists 
in computing the physical quantities in the "regular" regime as < e/b; the universal e-singular 
factors are then removed into non-perturbative quantities, and finally by analytic continuation 
the physical case at e = is recovered. 



2.4 Qualitative behaviour of the solutions 

Before embarking upon the resolution of the colhnear model equations (171[T2i) . it is instructive 
to estimate the qualitative behaviour of the solutions by using well-known methods [5] in the 
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context of high-energy QCD. Particularly important is the factorization property which allows 
one to split the unintegrated gluon density f{t,to) into a perturbative and a non-perturbative 
part, provided the "hard scale" t ^ > t\ is sufficiently large: 

f{t,to) = /pt(t)/np(to) X [1 + O (e-*)] , (15) 
up to terms exponentially suppressed in t (higher- twists). In turn, the perturbative factor 



/pt(t)~exp|y"*7(a(t')) dt' 



(16) 



is given in terms of the gluon anomalous dimension 7(a(t)) determined by the small-a; equation 

1 = a(t)x(7) , (17) 

where x is the eigenvalue function of the integral kernel in eq. ([I]). 

In this coUinear model, the eigenvalue function in eq. ([9]) provides two solution to eq. (ITTll 



7±(«) = \ , (18) 

the perturbative branch being the one with minus sign: 7- (a) = a + (9(a^). At large t, the 
running coupling saturates at the UV fixed point a(+oo) = 1/B = e/bu, so that the large-t 
behaviour of the unintegrated gluon density is given by 

/(t»to) ~$^c,(to)e*^^('^(+°°)) = c+(to)e^('+^) +c_(to)e^('"^) • (19) 
i=± 

According to the value of B we expect two kinds of asymptotic behaviour: 

• For B > 4 the square root is real and positive, the UV regular solution corresponds to 
the perturbative branch 7_ of the anomalous dimension and we must reject the (UV 
irregular) solution which diverges more rapidly: c+ = 0. 

• For B < 4 the two exponents are complex conjugate, and the gluon density becomes 
oscillatory at large t. It is not possible to distinguish an UV regular solution, and one 
has to determine the coefficients c± by analytic continuation in B from B > A. The fixed 
couphng {B = 0) solution belongs to this class. 

The above results will be also obtained in a more rigorous way in sec. 14.11 when treating the 
running-coupling equation. 



3 Collinear model with frozen coupling (b = 0) 

Since the properties of the solution of the collinear model and its connection with the solution 
method of ref. [1] are more easily illustrated in the fixed coupling case, I start considering the 
integral equation ([71) with 6 = 0. 
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3.1 Solution in momentum space 

The presence of the exponential factor in front of the r.h.s. of eqs. flUllTl) spoils scale invariance, 
therefore the determination of both eigenf unctions and eigenvalues of the integral operator by 
means of standard techniques is not possible. It turns out, however, that one can exactly 
solve eq. In fact, by differentiating it twice with respect to t, we obtain the second order 
differential equation (the e-dependence of / is understood in this section) 

/" - (1 + 2e)f' + [eil + e) + ae^/ = -ae^*« 5(1 - to) , (20) 

which can be recast in a more familiar form if we introduce the variables 

^ = i , ^ = 2r]V^^ , fit, to) = z^+^-J{z, zo) , (21) 



thus obtaining 

+ z-J' + [z^ - rf]'J = -^^^^ . (22) 

2r]z^ 



In the l.h.s. of eq. ( l22i) one recognizes the differential operator defining the Bessel functions 
J±r]{z) and Yrj{z) as solutions of the corresponding homogeneous equation. 
The general solution of eq. (122|) has the form 

3^{z, Zo) = ci{zo)3'i{z)e{zo -z) + cu{zo)MzMz - ^o) (23) 

where 5"/ and 3^u denote respectively the IR-regular and the UV-regular solutions of the homo- 
geneous equation, while cj and cu are zo-dependent coefficients to be determined by the two 
conditions of continuity of 5" and discontinuity of 8^3^ at z = Zq: 

\im 3'{z,zo)-lim3'{z,zo) = cu{zo)3'u{zo) - cj{zo)3'i{zo) = (24a) 

1 

2^ 



lim dz3'{z,zo) - lim dz3'{z,zo) = cu{zo)3''jj{zo) - ci{zo)3''j{zo) = - - — = -N{zo) . 



(24b) 



By solving the above linear system one obtains 



Zo) = [3^jiz)9^u{zo)Qizo -z)+ 3^u{z)3'iizo)e{z - zo)] , (25) 

where W = 3^u3^'i — 3^'u3^i is the Wronskian of the two solutions of the homogeneous equation. 
It remains to determine 5"/ and 3^u, each being a linear combinations of, say, and Y^: 

^s{z) = c^PUz) + cf\{z) , {s = I, U) (26) 

(the absolute normalization is irrelevant). From the asymptotic relations 

Xriiz) z"^ {z^ 0) Jriiz) ~ z-^^"^ cos{z + 0i) (z ^ +oo) (27a) 

Y^{z) ~ z-'^ {z 0) Yr,{z) ~ z-^^^ cos{z + 02) {z +oo) (27b) 

it is clear that the IR-regular solution is 5"/ oc J^, since it vanishes more rapidly than any 
linear combination containing Yrj{z) when z ^ with rj > 0. On the other hand, the UV- 
regular solution cannot be determined in this case of 6 = 0, because of the identical asymptotic 
behaviour (up to normalization and phase) for z +oo of all solutions in eq. ( l26l) . However, 
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the UV-regular solution can be unambiguously determined in the formulation with running 
coupling (cf. sec. 14.31) . and in the 6^0 limit it reduces to 3^u{^) = '^ni^)- conclusion 

3^j{z) = J^{z) (28) 
= Y^{z) = [cosiTT7])Jr,iz) - J-r,iz)]/ sm{nr]) (29) 
Wiz) = Y,J'^ - J^Y; = -2/7TZ , (30) 

whence 

S'iz, zo) = --^ [J^{z)Y^{zo)e{zo -z)+ Y^{z)J^{zo)e{z - zo)] . (31) 
^rjZQ 

It is possible to show that 3^{z, zq) in the previous equation obeys also the integral equation 
3.2 On-shell limit and perturbative expansion 

It is important at this point to check the explicit solution in eq. (1311) with known results of the 
literature. The perturbative expression for the GGF S{k, fco) in dimensional regularization was 
given in [7] (3.3) for an on-shell {ko = 0) initial gluon. In terms of the dimensionless density fe 
their result reads 



feit) = ae 



et 



m=l j=l 



(32) 



for a generic integral kernel with eigenvalue function x(7,e). 

The on-shell limit fco <^==^ to ~oo of the unintegrated gluon density f^it, to) at fixed 
t, e, ds is finite, and can be obtained from eqs. fl21ll3ip by exploiting the asymptotic behaviour 
of Bessel functions Jr^^z^) ^ (2;o/2)''/r(l + rj) for zq 0, whence 

m . ,_nn,^/.(M„) = .'^^Ito -^M^o)Y,M = -^^^ (i)'"i',(^) . (33) 

In words, the on-shell unintegrated gluon density is equal to the UV regular solution of the 
homogeneous differential equation with a proper normalization. 

In order to compare the solution flH^ with the perturbative expression one has to 

expand the r.h.s. of eq. (l33l) in series of z'^ ~ a. By rewriting as a combination of 
according to eq. (l29ll . and then using the ascending series [11](9.1.10) 

^M-(i)^f: J;;^:L) - 

m=0 ^ ^ 

one obtains 

fsit) = ^r(r, + r)sin(vrr,) ilY^^ ^^-^^'^ ' '^^^^^^ ^^^'^^^ (35) 

- (-r/W*)'" r(l-r^) _ . .^^^ . 2, . j-y^ae^T 
^ m\ r(l -r^ + m) 2^ {I + + m) ' 

m=0 ^ I ' J „=0 \ ' I ' I 

The first term in the r.h.s. of eq. exactly reproduces the perturbative result ([S2D, since for 
m> 1 



rn 1/2 1 ''^ 
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The second term of eq. (I35l) provides contributions of order (ae'^*)''^^^'" = a^+"^+'^/'^e[^+i'^+"'^)^]^ ^ 
each being outside the domain of the kernel and therefore out of the reach of the iterative 
procedure. Furthermore, this term is strongly suppressed ~ when e —>■ with respect to the 
perturbative one. Therefore, it is possible to correctly compute the perturbative coefficients to 
any order m provided e is sufficiently small {e < 1/m). In the limit e — > the perturbative 
solution agrees with the exact one to all ordersH 

As final remark, the series in eqs. fl32ll35p converge for all t G M, as one can check from the 
7 — > +00 asymptotic behaviour of xil) ~ 1/7- 



3.3 Integrated gluon densities 

The major issue this paper is devoted to, concerns the MS ^ Qq scheme- change, namely the 
relation between gluon densities and anomalous dimensions in the two factorization schemes. 
In the coUinear model, the off-shell integrated gluon density defined by 

g,{t,to) = I d^+^-fc' uj9{k\ko)e{e - k'^) = 1 + r dt' /,(t',to) (37) 

J J — oo 

can be computed in closed form (app. lA.ll) . and for t > to reads 

9e{t,to) = -7r^(^-^^ Jr,izo)Yr,+iiz) (t > to) . (38) 

Note the remarkable fact that g, like /, is factorized in the t- and to-dependence. 

The Qo-scheme gluon is given by the e ^ limit of the above expression, yielding (app. IA.2I) 

(t, to) = lim «,(t,to) = ^=^— -exp [7(a) (t - to)] (39) 
£^0 y/i — 4a7(a) 

whence one immediately derives the Qo-scheme anomalous dimension (a dot means t-derivative) 

7(«°)(«)= lim \q^;"\ =i{o) (g^dtg). (40) 

It is interesting to note that the on-shell limit of the integrated gluon density, i.e., the gluon 
density in dimensional regularization 

71 f Z \ ^"'"^ 

*W"3"oo*<'-'»' = -r(^TT)(2) ''"W 

provides the same effective anomalous dimension (app. IA.2p 

7efT(t) = hm ^ = hm ^ = hm — = 7(a) , (42) 

which means that the two limiting operations to —>■ —00 and e ^ commute. Actually, in this 
model this is a trivial consequence of the factorized structure of the gluon density (|38l) in its t 
and to dependence, which causes the ratios g/g in eqs. ( HOll and (H2l) to be to-independent. 



^These conclusions are valid in the ofF-shell case (^o G K) too, but for sake of simplicity they have been 
presented only in the on-shell case. 
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The relation with the MS-scheme anomalous dimension is obtained as follows. From the 
£ — > asymptotic behaviour of the on-shell gluon density (app. IA.2p 



one identifies the exponential in eq. (H3l) as the MS gluon density ^''^^^•'(t) jf] since it sums all 
and only e-singular terms up to the scale = /i^e^*. The MS anomalous dimension is then 
computed from the logarithmic derivative 

•(MS). X 
Qe 'it) 

and coincides, in this case of 6 = 0, with the Qo-scheme anomalous dimension, in agreement 
with refs. [7] and [1]. 

The coefficient function i?, on the other hand, is finite in the e — *• limit, and is given by 
the product [1] R = !N"9l, where 

-i 

N(a) = , = ^T7-^^-^^ (46) 

7(a)\/=?(^ 7(a)|l-4a]V. 

is the fiuctuation factor of the saddle-point estimate introduced in [1] (cf. also sec. 13.41) . while 
9^(a)=exp|^'^"^^d7| , x{l,e) = Xoil) + exiil) + O {e') (46) 

originates from the e-dependence of the eigenvalue function. Since in this model x is indepen- 
dent of £, Xi = 0, = 1 and therefore R = 'N, in agreement with eq. fHHj) . 

3.4 Solution in 7 space 

Having the solution of the integral equation at our disposal, we are ready to check the validity 
of the procedure suggested in ref. [1], at least in this simplified model. I start reviewing the 
main steps of that procedure. 

1) We introduce for the unintegrated gluon density /e(t,to) an integral representation of 
Mellin-type: 

fe{t,to)= I ^e^7e(7,to). (47) 

2) In 7-space, the integral equation ([7]) is thus recast into the finite difference equation 

A(7 + e, to) = ax(7)e-^*° + ax{l)fe{l, to) • (48) 

3) The finite difference equation (HHl) is solved in terms of a Laurent series in e, so as to 
provide the following expression (cf. [1], sec. 2 and eqs. (C.1,C.2)) for the on-shell unintegrated 
gluon density: 

m = Q{a, 5) ^ ^ e^* exp | ^ L{j') dj' - + g .-11^-^) (7) | , (49) 



■^Due to the particular definition of as in eq. ^ which includes e-dependent factors, eq. defines a 
"modified" minimal subtraction scheme, related to the customary MS and MS schemes by a finite scheme 
change. These details are unimportant for the purpose of this paper. 



11 



0.5 




Figure 3: Singularity structure of the Mellin transform fe{'y,tQ) in the complex '-/-plane. The 
shadowed region corresponds to the convergence strip of the Mellin transform; The crosses 
indicate the position of the singularities; the circles show the location of the poles of the terms 
in also shown are the original integration path C in eq. and the deformed contour C 
used in eq. (STi) . 



where Q = a/ a/2'!Te[l + O (e)] is a normalization factor, -^(7) = log [ax{"y)), -^^ = d!^L, and 
the coefficients i?„ denote Bernoulli numbers. 

4) The solution is determined by assuming the existence of a saddle point on the real axis, 
whose steepest descent direction lies on the real axis. 

Let us now analyze each point in turn, in the context of the coUinear model. 
1) Concerning the existence of a Mellin representation for the solution / of the integral 
equation ([71), the asymptotics in eq. (1271) guarantee that the Mellin transform / is defined 
in the strip 1/2 + 35/4 < 3^7 < 1 + e for all £ > 0. Explicitly, / is given in terms of 
iF2(m; u+ 1,v] —Zq/4:) sums, as follows: 

y + jdte-'r'Mt,to)^ft\l,to) + fl+\j,to) (50a) 

00 / 2 I A\k 

/i^H7,to) = vr (1)"'^ J,(.o) {-{a^'rcot{7.ill) ^^^ + (50c) 



sin(7rr7) ^ (^ - 1)! 



1 {zl/AYcoi^{TTr]) 



{k-r]-f)T{k-T]) [k + 7]{l--f)]T{k + r]) 



I will show now that, with a proper choice of the contour C, only the ffist term of f^^^ in 
eq. (I50cp contributes to the inverse Mellin transform fl47l) for t > to — the relevant region for the 
on-shell limit. Notice that the analytic continuation of / defines a meromorphic function whose 
singularities are just the simple poles of f^~^ at 7 = 1 + fee : = 1, 2, ■ ■ ■ , as shown in fig. 13 
Actually, f^^^ is holomorphic in the whole plane 7 G C, since the poles at 7 = fee : A; = 1, 2, ■ ■ ■ 
stemming from the ratio cot(7r?77)/r(?77) in the ffist line of eq. fl50cp are exactly canceled by 
those in the sum on the second line; furthermore, the poles at 7 = 1 + ke : k = 1, 2, ■ ■ ■ 
stemming from r(l + ?7(l — 7)) in the first line are also canceled by those in the sum on the 
second line. 

It is convenient to compute the inverse Mellin transform separately for the (— ) and (+) 
pieces. In the integral of e^^ f^~\'j,to) one can close the contour path to the left {t — to > 
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0), without crossing any singularity, thus obtaining a vanishing contribution, as expected. 
Considering now the integral of e'^*/^"^-'(7, to); one is not allowed to close the contour either to 
the left or to the right, because the factor e'^''*~*°'* in front of the sum grows for 9fJ(7) — »■ +00, 
while the ratio of gamma-functions in the first term grows with I7I for 3^(7) < 1/2 + e/2. 
However, by folding the contour C ^ C so as to let it cross the real axis at some value 70 < £ 
(remember that f^^^ has no singularity), and then computing the two contributions of eq. f lSOcI) 
separately, one obtains a vanishing integral from the second line, because the contour can be 
closed to the left without crossing any singularities. 

To summarize, with an integration contour C crossing the real axis at 70 < £ and going to 
infinity with 3^(7) > 1/2 + e/2 as in fig. [3], only the first term in eq. flSOcI) contributes in the 
7- representation (H71) for t > to- 

By performing the on-shell limit I end up with 



Mt)= hm /,(t,to) 



/ ^ e-^(a,^)-cot(vr,7)^^^^^#T^ - / ^ e^'fsh) , (51) 



r(l + 7]) Jc> 27ri 

which is just the Mellin-Barnes representation [11] (9. 1.26) of the Bessel function in eq. 

2) It is straightforward to check that the on-shell Mellin transform fe{l) in eq. (|5T1) obeys 
the homogeneous difference equation 

A(7 + £) = ax(7)/.(7) (52) 



analogous to eq. [1](2.11). With some more effort, one can show that the off-shell expression (1501) 
obeys the inhomogeneous difference equation (HHl) . 

3) The third issue concerns the validity of eq. fH9|) . By explicitly computing the integral 
and the derivatives of L{'-y) in the coUinear model 

I L(7')d7' = 7 log(a) - 7 log(7) + (1 - 7) log(l - 7) + 27 

L(-)(7) = (m - l)![(-l)-7-"^ + (1 - 7)-"^] (m > 1) , (53) 
the exponent within curly brackets in eq. ( H9l) becomes (-8201+1 = : m > 1) 
^(7) = r7[7 log(a) + 27-7 log(7) + (1 - 7) log(l - 7)] (54) 

- \ log(a) + \ iog(7) + \ iog(i - 7) + E wirTT) i - - 1^^]""} • 

m=l ^ ' 

The sum in the above equation is typical of the asymptotic expansion of the logarithm of the 
gamma- function [11] (6. 1.40). In fact, by comparing eq. (15^ with the asymptotic expansion 

logr(?7(l - 7)) - logr(r77) ^ r/ [(1 - 27) log(r/) - 1 + 27 - 7log(7) + (1 - 7) log(l - 7)] 

+ 2 log(7) - 2 log(l - 7) + E 2m{2m- I) ^ '^^^ " ^^l"" ' t^^]""^ ^^^^ 

m=l ^ ' 

one gets 



r(i + - 7)) 
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Im(y) 



Im(r) 




Figure 4: a) The imaginary part of the integrand in eq. ( f3I]) showing the singularities on the 
positive real semi-axis; in yellow a sketch of the fastest convergence path, b) asymptotic limit 
of the integrand showing the discontinuity ( f5^) on the real axis with a peak around the saddle 
point value ( f59l) . 

Apart from the irrelevant normalization factor ^ j ^farf^^^ eq. fl56l) agrees with the integrand 
in eq. flSip . when one takes into account that the ri — * +00 asymptotic expansion in powers of 
l/?7 of cot(7r?77) is a numeric constant (=Fi according to the sign of ^7). 

4) The last step is to evaluate the integral in eq. ( 15T1) in the large-?] limit. It turns out that, 
for small values of e and values of ae*^* < l/x(l/2) = 1/4, the fastest convergence contour path 
surrounds the interval < 3fJ(7) < 1/2 (cf. fig. H]) at a distance decreasing with e. The main 
contribution to the integral comes just from this region (parts B and D in fig. S]). In the limit of 
vanishing e, the string of poles at 7 = accumulates into a branch-cut at 7 g]0, +cx)[. In fact, 
while the ratio of gamma-functions is regular at 7 > also in the r\ -|-cxd limit, the cotangent 
cot(7rry7) — isign(SJ7) becomes discontinuous across the real axis with a jump equal to — 2i. 

Therefore, neglecting the contributions to the integral in eq. flSlI) from the parts A, C and 
E of the contour path, the contributions of B and D amount to the integral in 7 g]0, l/2[ of the 
discontinuity of the integrand, which can be easily obtained by replacing cot(7r?77) with — 2i. 
One obtains 

where use have been made of eq. 

Some remarks are in order. Firstly, by expanding in e — ^ the prefactor in the r.h.s. of 
eq. (1571) ^er'^rp^'^ jVij] + 1) = \/ a/27re[l + O {e)\ eq. (H9l) is correctly reproduced. Secondly, 
the integrand, being a discontinuity of a solution of the difference equation (|52i) . is itself a 
solution of the same equation. Thirdly, the integral representation (1571) of the on-shell density 
/ uses an integration path lying on the real axis. 

The last remark explains the possibility of having a stable saddle point in the real direction, 
despite the fact that our original integral in eq. ( |471) involves a real analytic integrand and an 
integration contour C parallel to the imaginary axis. In fact, according to the analysis of [1], in 
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the 1] —>■ oo limit the leading part of the exponent in eq. (1571) is given by 

7t + 5(7) ~ 7t + r/ r L{j') dy + O {rP) , (58) 
Jo 

having considered t a possibly large parameter. The saddle point condition is (cf. eqs. ( IT71IT81) ) 



^t + logM7)] =0 ^ 7(ae-O=^-^^-^^ = ae-* + 0((ae-T) (59) 
_X'(7) < 2 

and is fulfilled when 4ae'^* < 1 so that < 7 < 1/2. The saddle point behaviour of the 
discontinuity of the original integrand /(7) is apparent in fig. |H The final result is 

1 {IP 1 ae^^ { 1 f""^^^ da 

/.(t)-^^— exp|7t+-y^ log[ax(7')]d7'| = |^3^^exp|-y^ -7(a) 

(60) 

and exhibits the factorization of the 1/e collinear singularities, thus allowing us to derive the 
relation between the MS and Qo gluons, as explained in the previous section. 

In conclusion, the analysis of the frozen-coupling collinear model provides analytic expres- 
sions for the gluon densities and anomalous dimensions in both MS- and Qo-schemes which 
agree with the results of ref. [1], sec. 2. In particular, the explicit expression of the Mellin 
transform fei"^) offers a concrete test of the asymptotic series representation ( H9il . and also 
allows us to understand the relation between the original Mellin integral ( l5Ti) and the real-axis 
integral (1571) . the latter being the basic tool to prove (by saddle-point estimate) the factorization 
of collinear singularities 



4 Collinear model with running coupling {h > 0) 

In this section I shall extend the collinear model to the more realistic situation of running 
coupling. I shall show that most of the analysis preformed in sec. Ofor the fixed coupling case 
can be carried out with running coupling too, with analogous results. 

4.1 Asymptotic behaviour of the solutions 

Like the 6 = case, the solution of eq. ffT^ obeys a second order differential equation: 

(1 + ABe'')f" -{1 + 26 + ABe'')f + [e{l + e) + Ae^y = -yle"*°5(t - to) . (61) 

In order to characterize the IR and UV regular solutions of the homogeneous equation, I 
first determine their large-|t| behaviour. This can be accomplished by rewriting eq. ( I6T1) in 
Schrodinger-like form and then using the WKB approximation. In detail, by letting 

we obtain for h the Schrodinger equation 
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Figure 5: The potential V{t) of the Schrddinger equation ( f63]) in the three running- couphng 
regimes: Rl (sohd-red), R2 (dashed-green) and R3 (dotted-blue) . 



The WKB approximation to the solution of eq. fl63|) . written in terms of the wave- number 
K = VV, reads 



hit) 



-^exp{±/.(t)dt} 



(64) 



and yields, when inserted into eq. (lU^ . the two possible asymptotic behaviours of f{t). 
In the IR region [t — oo) we have 



exp 



1 1 . 
- + e± - ]t 

2 2 ' 



(65) 



and one identifies the IR regular solution as the one which vanishes more rapidly, i.e., // 



In the UV region [t — > +oo) we have 



A ~ e*/2 , V 



4 B ' 



exp 



(66) 



and the solutions can have an exponential or oscillatory behaviour according to whether B is 
greater than or less than 4. In the former case, one again identifies the UV regular solution as 

the one which vanishes more rapidly, i.e., fu ~ exp 



ill 



1-5 1« 



it is convenient to 



According to the value of 5 — and always considering a < 1/^.-. 
distinguish 3 regimes where the "potential" V is qualitatively different (cf. fig. [5]): 
Rl: 0<5<4 

The potential is regular for all values of t G M, since A > 0, and its UV limit is negative: 



V{+oo) 



1 1 



(67) 



As a consequence, in the UV region the wave-number k is pure imaginary and all solutions 
share the same oscillatory behaviour, up to a relative phase. The UV regular solution is thus 
undetermined. 

R2: A<B < 1/a 

The potential is regular for all values of t G M, since A > 0, and its UV limit V{+oo) is positive. 
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Hence, the wave-number k is real throughout the whole t range and the UV regular solution is 
uniquely determined, as discussed above. 
R3: 1/a < B 

The potential is singular at the Landau point of eq. dHl), since A < and the denominator 
of V vanishes at t = ^a- This singularity might prevent the existence of a global solution for the 
integral equation ( fT2l) . Nevertheless, The UV regular solution of the differential equation ( 16T|) 
can be unambiguously identified. 

In the next section I shall explicitly determine the solution of the integral equation (1121) 
in the intermediate regime R2, leaving to subsequent sections the analysis of the regime Rl 
relevant in the limiting case 6 — > and of the regime R3 where the physical situation £ ^ at 
fixed h is recovered. 



4.2 Solution in momentum space 

By introducing the new variables 

a{t)B 



c 



-A5e 



1 - a{t)B 



(6J 



the integral equation ([6]) becomes 

1 



-c 



-Co 



-C')"'-^(C',Co) dC 



^(C, Co) = ^^tzq I - + 0(C - Co 

By differentiating the above equation twice with respect to C yields the differential equation 

V 



(69) 



C(l-C)^"+[(l-r]) + (r]-3)C]^' 



B 



r/ + 1 



B 



S{C - Co) 



(70) 



whose homogeneous version is just the hypergeometric differential equation with parameters 
u, V] w given by 



M, V 



1-? I 1±a/1-- 



w 



1] = U + V — 1 . 



(71) 



Let me now consider the regime R2 in which 4 < B < 1/aso that u < v are both real, A > 
and C < is a decreasing function of t (cf. eq. ( |68l) ). At variance with the 6 = case, both IR 
(C — > 0^) and UV (C — > —oo) regular solutions of the homogeneous differential equation (!70l) . 
are unambiguously identified, as pointed out in sec. 14. 1[ Explicitly 



j-^(C) = c^(-C)-%Fi 



2 — u , V 
V — u + 1 



cu 



r(2-M)r(i-M) 

T{l-w)T{v-u + l) 



- u , 2 — V 
2-w 



w[j^u,j'i] = {w-i){-cr'"i^-cy 



(72) 

(73) 
(74) 



By repeating the steps outlined in sec. 13.11 the conditions of continuity of at C = Co and 
discontinuity of the first derivative — A^(Co) = '7/[-BCo(l — Co)] provide the solution of eqs. (!69l) 
and (CO]): 

•^(C, Co) = 7^7^ [^,(C)-^f/(Co)e(C - Co) + ^c/(C)-^/(Co)e(Co - 0] , (75) 
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The integrated gluon density defined in eq. (1371) can be computed in closed form, and for 
t > to reads (app. IA.3P 



9e,b(ty to 



Cu 



B{l-v){-Co)^ 



2-u, v-1 

V — u + 1 



(76) 



showing also in this case a factorized structure. It is interesting to note that the dependence 
of the equations and their solutions on t occurs only through (, and because of eq. flUHl) . only 
through the combination a{t)B. 

The on-shell limit Co — ^ 0^ of both integrated and unintegrated densities can be easily 
computed by noting that ^/(Co)/(— Co)'' 1, whence 



fe,b{t) 
9e,b{t) 



B 

Cu 



V 



B(l 



2-u, v-1 
V — u + 1 



(77) 
(78) 



The comparison of the unintegrated gluon density / with the perturbative solution can be 
obtained rewriting the UV regular solution JF^ as the sum of two hypergeometric functions 
with argument ( by means of the inversion formula [11] (15.3.7) 



^c/(C) = 2F1 



U , V 

w 



T{v)T{v- l)T{l-w 
and then using their series representation [11](15.1.1). The first term yields 



(79) 



B 



u , V 
w 



= Ae'' h +^(Ae"*)" JJ 

t n=l k=l 



u + k-l){v + k-l){-B) 



{w + k- l)k 
Ae^' \ 1 + ^(Ae^*)" n ixike) - B] \ , 



(80) 



n=l 



k=l 



and provides the perturbative expansion in terms of the parameter yl = a/(l — aB) and of the 
"effective" eigenvalue function x{l)~B relative to the kernel defined in eq. (|T3|) . The additional 
contribution to f(t) due to the second term in eq. fl7^ is of order {—Qj^j ~ (— C)^"*" ~ 
(ae^*)*^''"'"^-* and cannot find place in the iterative solution, because it does not belong to the 
domain of the kernel K (cf. sec. 13.21) . However, in the e — limit, the perturbative solution ( IHOi) 
agrees with the exact one ( 1771) to all orders. 

As last remark, the domain of convergence of the series in eq. (IHOj) is finite (e^* < = 
|1 — e/foasl), unlike the 6 = case. 



4.3 Fixed coupling limit h — ^ 

It is important to study the limit 6^0 because, as already mentioned in sec. 13.11 it is not 
possible to determine the UV regular function in the frozen coupling case. Actually, we saw 
in sec. 14.11 that this problem is also present at 6 > when B < A (regime Rl). Therefore, I 
shall first derive the expression of the UV regular solution J-'u for B < 4 and then compute its 
limiting result at -B = 0. 

Clearly, 5 = 4 is a point of non-analyticity for the coefficients u, v and therefore also for 
the functions J-'j, Tu- When B < 4, the coefficients u, v become complex conjugate. It is easily 
verified that the IR regular solution jFj remains real. On the other hand, the expression for 
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Figure 6: Unintegrated gluon densities (times e~*/^ j in various running- coupling regimes, with 
parameters a = 0.2, r] = 4.5 and = 0. On the left side: lowest perturbative order (dash- 
dotted-black) , fixed coupling (solid-red), regime Rl (dashed-purple) , regime R2 (dotted-blue); 
on the right side: again regime R2 (dotted-blue), boundary between R2 and R3 (dashed-green) , 
regime R3 (solid-cyan) and e — > limit at fixed b = 1.2 (dash-dotted-brown). The last two 
curves diverge at their Landau points. 



Tjj in eq. fl72l) yields two different (complex) results depending on the choice '^u = —^v either 
greater or less than zero — corresponding to an analytic continuation from 5 > 4 to i? < 4 in 
the complex variable Bq, from above [Bq = B + iO) or from below [Be = B — iO). 

Note that S = 4 is by no means a critical value for the coefficients of the differential 
equations ( l6T]l70i) : it only separates the regimes of positive and negative effective potential in 
the UV region. Since nothing prevents the existence of a real solution, it seems reasonable to 
define the UV regular solution at < 4 by taking the average of the two analytic continuations: 



J^uiC; B<4) = - [J^uiC; B + iO)+ MC; B - iO)] = 3?^t;(C; 5 ± iO) . 



(81) 



Of course, the definition flHTl) joins continuously with the original definition fl72|) at 5 = 4. An 
explicit expression of the UV regular solution for all i? > can be obtained by applying the 
prescription ( ISTl) to eq. ( I79l) : 



:Fu{C) = 2F1 



w 



u , 2 — V 



w 



c 



(82) 

The B ^ limit, at fixed a, e, t is then performed by exploiting the series representation for 
the hypergeometric functions and the Stirling approximation for the ensuing gamma-functions: 



2F1 



2F1 



u , V 



w 

2 -n, 2-v 
2-w 



C 



3fJ 



V{2-u)V{l-u) 
V{v)V{v-l) 



r(i + ^) {^''U'^) 

^-^'^]' = cos(7rr^) 



(83a) 
(83b) 
(83c) 
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Note in particular the cosine term in eq. f l83cl) stemming from the real part of the complex 
exponential: it is exactly the "relative weight" between and needed to build up the 
Bessel function of the second kind F^, according to eq. fl^ . In conclusion, the substitution 
of the expressions ( l83l) into eq. f TfSj) yields the fixed-coupling result ( !3T|) . when the proper 
normalization factors between 3^{z) and J-'iC) are taken into account. 

4.4 Integrated gluon densities 

In this section I show how to derive explicit expressions for the integrated gluon densities 
and anomalous dimensions in the MS-scheme and Qo-scheme in the running coupling case. In 
particular, in this model it is confirmed the claim of ref. [1] that the running coupling corrections 
to the MS anomalous dimension 7*^^^) (a(t), 6) are provided to all orders by the e-dependence 
of the eigenvalue function x(7,e), according to the relation 

1 = a{t)x{l^'''\ a{t)bu;) , (84) 

where e has been replaced by a{t)huj. 

4.4.1 MS-scheme 

The MS gluon for 6 7^ is defined in two steps. First, at fixed h and e, one perturbatively 
{ots 0) computes the integrated gluon in dimensional regularization. This implies that the 
calculation is naturally performed with ots < s/h^ i.e., a < 1/B, which corresponds to the 
regimes Rl and R2. Then, all ensuing IR singularities appearing as poles at e = are isolated 
and factorized into an IR-singular "transition function", to be identified with the MS gluon 
density g^^^\t). 

The important point characterizing the MS-scheme is the factorization of e-poles in the 
form 

^(^s)(t) = exp||_* dr7(^^)(a(r),6)}=exp|i£^*^-^^^7(^'Ha,6)| , (85) 
where the MS anomalous dimension function 

^(MS) ig 

required to be e-independent. The first 
integral in eq. flH^ . which is singular for e — > because of its IR lower bound, defines the 
MS anomalous dimension. The second integral is obtained by changing integration variable 
according to eq. ( |TOl) . and is more suitable for comparison with perturbative calculations. If 
eq. (!85ll contains all IR singularities, the integrated gluon density (178|l of the coUinear model 
can be decomposed in the product 

f 1 dfl 1 

9eAt) = Re{a{t), b) exp I - ^(Y^^ T^'^'Ha, [ , (§6) 

where the coefficient function R is regular at e = 0. 

The above expression suggests a method for extracting the MS anomalous dimension. One 
observes that the integrand in the exponent is singular at i? = 1/a. On the other hand, for 
B — >• l/a(t), no singularity occurs in the off-shell functions /^^^(t, to) and ge,hiti to)- This signals 
that such singularity in the on-shell limit is connected with the infinite evolution of r from t 
to to = —00, and therefore it affects only the exponential factor, while no such singularity is 
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expected in the coefficient function R. Therefore, if we take the logarithmic derivative of g with 
respect to a{t) and subsequently the limit B —>■ l/a{t) from below, we obtain [at = a{t)): 



hm [l — atB)— 

B^l/at oat 



lim (1 — afB) 



da,R ^ ^^^^\auBe/u:) 
R Sdti^ — citB) 



eat 



Since the limit can be computed at any a{t) and e, the above formula enables us to deduce the 
full functional dependence of both a and b. In this model, from eq. flTHj) we get 



dlogg dC dlogg 



-B 



dat 



dat d( 



;i - atBY 



1 -V 



1 {2-u){v-l) 



2F1 



U , V 



V - u + 2 



C C v-u+l 



2F1 



2 — M , V — 
V — u + 1 



where eq. ( l68ll and the differentiation formula for hypergeometric functions [11](15.2.1) have 
been used. By performing the limit B 1/at, which implies l/( ~ atB — 1^0, the second 
term within square brackets does not contribute, and we end up with the simple expression 
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(MS), 



at,b) = eat hm (1 — atB) 



dhgg 



dat 



ed-v) 



\ B=l / at ,e=atbu) 



(89) 



e=atbuj 



which coincides with its fixed coupling (6 = 0) counterpart computed in eqs. (jHj) and (l59l) . This 
is a non-trivial result, since it shows that the e- independent kernel (|H]) provides a 6-independent 
MS anomalous dimension, according to eq. (15^ . 

Actually, it is not difficult to extend the coUinear model to e-dependent kernels and check 
eq. (!84l) in situations where 

^(MS) ig 

explicitly 6-dependent. For instance, by considering an 

e-dependent kernel 



K^"\r,e) = H(.)[0(-r)e«(^)^ + Q{r)] , x{l,e) = E{e) (^^ + ^(^) 



(90) 



where S and are regular functions of e, one obtains the same type of differential equation 
and hypergeometric solutions. Analogous expressions hold for the gluon density, with the 
replacements A —>■ AS, B —>■ B /S and with new parameters 



M,'y = l- y|^l±yi-^| , w = l-^r] . 
The MS anomalous dimension is then straightforwardly obtained (cf. eq. 



(91) 



(MS)/ M n M Cf h ^ 1 - - 4at E{atbuj) /^{atbuj) 
Y '{at,b) = e{l-v)\B=i/at;e=atbL, = ^{atbuj) . (92) 



On the other hand, the solution 'y{at,e) of the implicit equation 1 = atxil,^) (satisfying the 
perturbative condition 7(0, e) = 0) is given by 



-/ N n1 - yl -4atS(£)/^(e) 

l{at,£)=^{£) = e{l - v)\B=i/at (93) 



and exactly reproduces the anomalous dimension in eq. ( I92l) when substituting e atbuj 
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4.4.2 Qo-scheme 



The Qo-scheme gluon is defined by the e ^ hmit at fixed b of the off-shell 6-dependent 
integrated density ( I76l) . analogously to the frozen-coupling definition in eq. ( l39l) . In this limit, 
B +00 and we enter the regime R3, in which the variable ( is positive. Therefore, we 
must extend the expression ( 1761) of the integrated density from R2 to R3. At the separation 
point B = 1/a, however, the running coupling assumes a constant value a{t) = a, whereas 
the parameter A (eq. fill I) ) and the variable ( diverge. In order to avoid the singular point 
B = 1/a, we can change regime by means of an analytic continuation to complex B of the 
equations fl7m78p obtained in R2. According to whether R3 is reached from the upper or lower 
half of the Be = -B±iO complex plane, (— C)^ ~^ (^Pg±i7rp acquires a phase of different sign. This 
translates into a discontinuity of the analytic continuation of g at values of B > 1/a. For t > 
we obtain 



B{l-v){u- 

giito) = (1 - Co) 



2-u, v-1 
V — u + 1 



u , 2 — V 
u — V + 1 



1 

Co 



^ ) giito) , {t>to, B> 1/a) (94) 



M-v) r{2-u)r{i-u)r{u-v + i) 

T{2 - v)T{l - v)T{v - u + ly^ 



2 — u , V 
V — u + 1 



1 

Co 



(95) 



From this equation we learn that: 

• the analytic continuation of the gluon density is still factorized in its t and to dependence 
with the same UV (t-dependent) factor as in R2; 

• the discontinuity affects only the IR factor gi{to), because of the phase e='=^'^("~") in the 
second term; 

• the Landau pole shows up in the well known branch point of the hypergeometric function 
at C = e^(*~*A) = 1^ in both UV and IR parts. 

Our main goal, though, is to obtain the anomalous dimension, which is known to be inde- 
pendent from the IR properties of the theory, provided the "hard scale" t is large enough. In 
other words, the effective anomalous dimension 



7eff(t,io 



gs,b{t,to) 



f2-u,v 
1-v ^\v-u + l 






('2-u, V- 1 

y V — U +1 


9 



{g = dtg) 



(96) 



is expected to depend only on a{t) for t ^ tQ,tx, where eventually represents a cutoff 
that regularizes the Landau pole and gives mathematical meaning to the gluon density g. In 
this model, thanks to the factorization property of g, the effective anomalous dimension is 
independent of t^, hence independent on the details of the analytic continuations to R3, and 
needs not a regulator of the Landau pole. As a result, the Qo-scheme anomalous dimension 



7 



(Qo) 



{ait)) 



y y gb{t,to) 

hm hm — 

io^-oo £^0 gi^i^t, to) 



lim7efT(t) 



(97) 



is just the £ — > limit of the effective anomalous dimension. 
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In order to compute I rewrite the hypergeometric functions in eq. ( I96ll as [11] (15.3.5) 



2F1 



2 — u , V — 1 + n 
V — u + 1 



1 



l—n—v 



, V -1 , V - 1 + n X \ , - 
v-u + 1 wj ' (- = 0'1)- 

(98) 

In this form, only the third parameter and the argument of the hypergeometric function diverge 
in the 77 — > 00 hmit. Finally, by using the limit representation [12] of the Tricomi confluent 
hypergeometric function U 



1-C 



lim 2F1 



a , b 



z^t/(a, a-b + l,z) 



(99) 



and the four- dimensional running coupling a{t) = a/(l + abut) = l/buT, I obtain 



7(««)(a(t),6) =a{t) 



U{-l/buj,0,T) 
U{-l/buj,l,T) 



T = t + 



abuo 



(100) 



which is indeed a function of a{t) and buj only. Let me stress that the expression above resums 
the running coupling corrections of the Qo-scheme anomalous dimension to all orders in b. 

The Qo-anomalous dimension f llOOl) can be directly obtained in D = 4 dimensions as well. In 
fact, the differential equation ( l6Ti) at e = 0, in the T variable reduces to the Kummer's equation, 
whose independent solutions are the confluent hypergeometric functions iFi(l — l/6ti;, 2, T) and 
U{1 — l/buj,2,T) = U{—l/bu!,0,T)/T. The UV asymptotic behaviour identifies U(T) as the 
UV regular solution, i.e., as the unintegrated gluon density f = g — up to an IR-dependent 
normalization constant. In turn, U{1 — 1 /boo, 2, T) is the derivative of bu!U{—l /bu, 1, T), which 
represent therefore the integrated gluon g — up to the same IR-dependent normalization. 
Finally, the ratio g/g yields the Qo-anomalous dimension as in eq. f llOOp . 

An important check comes from the 7-representation [13] 



d7 
27ri 



d7 
27ri 



7) 



1/buj 



^(7) = / X(7) d7 = log ■ 



7 



-7 
(101) 

yielding, in the coUinear model, the confluent hypergeometric function U{l — l/buj, 2, T) — apart 
from a t-independent factor — as first observed by M. Taiuti in her degree thesis [14]. From the 
representation f llOip . by means of the saddle-point method, one can obtain the running coupling 
corrections at any given order in bu. In the 6 — >■ limit, the Qq anomalous dimension fllOOl) 
reduces to the frozen coupling value 7(a(t)) of eq. (HUj) . and coincides with the MS anomalous 
dimension 0891) . Starting from 0{b), the two factorization schemes provide different results. 
In particular, this collinear model predicts a 6-independent MS-scheme anomalous dimension, 
whereas the Qo-scheme contains non-vanishing corrections, which agree with the 6-expansion 
in D = 4 dimensions. 



4.5 Solution in 7 space 

In this section I show that the 7-representation (1471) is valid also at 6 7^ and similar conclusions 
as for the 6 = case can be drawn in the asymptotic small-e expansion of the gluon density. 

The Mellin transform fe,h{l^ to) of the unintegrated gluon density, defined in eq. fl50ap . exists 
in the strip 1/2 — 3?a/1/4 — 1/B < 3?7 < 1 + £ for all values of e and B. However, the structure 
of the singularities of /(7) is different from that in eq. ( l50l) . since now has two infinite 
series of poles at 7 = (1 ± a/1 — 4/i?)/2 — en : n G N. In particular, when B < A (regime Rl), 
these poles are located off the real shown in fig. d 
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Figure 7: Singularity structure of the Mellin transform fe{l,to) in the complex 'j-plane. Sym- 
bols as in fig. [3l 



Let me focus on the gluon density for t > to- By following the same steps as in sec. 
the function f^^\'j) can be decomposed in the sum of a product of elementary and gamma 
functions plus hypergeometric functions of type 3F2; by then deforming the integration contour 
C into C so as to cross the real axis at 70 < e (while leaving the new complex poles to the left), 
the 3F2 functions do not contribute to fe^lt > to). In the on-shell limit to ~oo we are left 
with 

r(l + r7(l -7)) 
X f^— -^[-cot(7rr/7) + 3fJcot(7rM)] , (102) 
r(777) 

where the real part in the last term descends from the real part used in eq. ( l82l) . With some 
trigonometric identities it is not difficult to prove that, for B < A, 

cr\ / \ -^r/\ /\T Sin(7r'?^) rj^oo \ r~A 

iKcotynu) = -[cot^TTM) + co^Tcv)] = . ~ exp —T^f]\ B ~ 

cos{-KT]) — cosh (vrr^w -| — l) 

(103) 

When computing the inverse Mellin transform fll02p along the deformed path C in the large-?] 
(fixed B) limit, the term with — cot(7r?77) —>■ ±isign(537) becomes discontinuous on the real 
axis and its main contribution is provided by the integral of such discontinuity in the real 
interval 7 g]0, 1/2[. On the other hand, the term proportional to ^cot{nu) does not develop 
any discontinuity, hence does not contribute in the parts B-C-D of the contour; furthermore 
in the remaining parts A and E it is exponentially suppressed with respect to the other term 
cx cot (71777), as indicated in eq. f llOSp . and therefore can be completely neglected. In conclusion 

^•■'") " grwrwre.) /*'^'^"'-''"'" r(« -i + .nOr(«-i + ,„)£ii±^^ , (io4) 

The main difference between eq. ( 11041) and its 6 = counterpart eq. ([57j) is the presence of 
the two additional u, f-dependent gamma-functions. The latter modify the analytic structure of 
the Mellin transform f{'~f) away from the real axis, but do not affect the mechanism generating 
the discontinuity in the £ — > limit. It is an easy exercise to check that in the b ^ case the 
integrand of eq. (11021) (without e'''*) and its discontinuity (I104p obey the homogeneous difference 
equation 

/,,,(7 + s) = Ax(7,i?)/.,6(7). (105) 
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Finally, in the regime Rl we can verify the validity of the Laurent- series representation (H9l) 
also for the running coupling case. In fact, by using in eq. fll04p the asymptotic expansion of the 
gamma-functions in terms of Bernoulli numbers, after some calculation one indeed reproduces 
eq. dlH]) with 

|7'7+l/2g-r? T(n/B) 

"=Vi^ r(»,+ l) (,/flA4-.,/B - i(7.B)=log|^xh,B)|. (106) 

In the 7] — > +00 limit, \/ A/2Tie and 

- Tfc / {; ^ " } 

agrees with eq. [1](4.1). 

To conclude this section, I have shown that the analysis of the collinear model can be 
carried out explicitly in the presence of running coupling, and provides analytic results for the 
anomalous dimensions in both MS- and Qo-scheme, which agree with the general results of the 
literature and in particular with the relations provided by sees. 3 and 4 of ref. [1]. 

5 Conclusions 

In this article I have considered a simplified version of the integral equation that determines the 
gluon Green's function in high-energy QCD in arbitrary space-time dimensions D = 4 + 2^. The 
kernel of the integral equation agrees with the true leading-log x BFKL kernel in the collinear 
limit, where the transverse momenta of the gluons are strongly ordered. This model has no 
phenomenological ambition, but embodies most of the qualitative features of the real theory, 
e.g., the kinematical symmetry in the gluon exchange, the leading- twist behaviour of the gluon 
density, the pattern of IR singularities and the running coupling. It is therefore a useful tool 
to check and better understand general results of the QCD literature. In fact, this model was 
already considered in Z) = 4 dimensions [9] for clarifying the transition mechanism between the 
perturbative, non-Regge regime and the strong coupling Pomeron behaviour. 

In the present formulation I have explicitly determined the gluon densities and their anoma- 
lous dimensions in two different factorization schemes: the MS-scheme, based on dimensional 
regularization, and the Qo-scheme, based on an initial off-shell gluon. The main motivation for 
this analysis stems from a previous work by M.Ciafaloni and myself [1] where we introduced 
a new method for solving the off-dimensional BFKL equation and for performing the mini- 
mal subtraction of the collinear singularities. The rather formal expressions we obtained and 
some sensible but unproven assumptions we made, could be strongly supported by an explicit 
non-trivial example where they are shown to be valid. 

This analysis attains this object. In fact, in 4 + 25 dimensions, the master integral equation 
is solvable in terms of Bessel functions (with frozen coupling) and hypergeometric functions (in 
the running coupling case). The results obtained here are then often compared with series and 
integral representations of ref. [1], showing their correctness and their domain of validity. In 
particular, it is clarified the mechanism by which the integral representation of the solution of 
the master equation — a real analytic function of the anomalous dimension variable 7 integrated 
along a contour parallel to the imaginary axis — is evaluated by a saddle point integral along 
the real axis. I also show that the iterative/perturbative solution to the master equation agrees 
with the exact one, at least up to the order n < higher orders n> 1/e of the perturbative 
expansion do not belong anymore to the domain of the kernel. Among the most important 



25 



results there is the confirmation of the formula (l8^ determining the MS anomalous dimension 
with running coupling from the e-dependence of the kernel, which in general was proved only 
up to O (5^) corrections [1]. 

On the whole, this model represents a useful tool for studying the mathematical properties 
and the qualitative features of the off-dimensional BFKL equation, even with running coupling. 
It supports the validity of the procedure [1] for determining anomalous dimensions in sublead- 
ing approximation, and encourages its application for extracting the leading-twist anomalous 
dimension at full NLx level. 
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A Integrated gluon density 

In this appendix I compute the off-shell integrated gluon density defined in eq. (1H7|) . both at 
fixed coupling and with running coupling. 

A.l Integrated gluon with frozen coupling 

At fixed coupling, it is convenient to use the z- variable introduced in eq. fl^Tj) : 

g,(t, to) = l + 2r] [ zo) dx , (108) 

where the unintegrated density 5" is given in eq. (l3Ti) . For t < to we have [11] (11.3.20) 

g^it < to) = l + 2r/ ~^/;lf°^ fx^^'Ux) = 1 - ^ (-\ J,+^{z)Y,{zo) . (109) 



For t > tQ, the integral in eq. fllOSp is conveniently split into 2 pieces 



27] [ x''+^3'{x,Zo) dx 
Jo 



Y^izo) x''+V,(x)dx + J,(zo) / x'^+%ix)dx 



TT 



2z^, 



(110) 

[^0^+1 {y,(zo) Vi(^o) - J,{zo)Y,+i{zo)} + J,{zo)z'^+X+,{z)] . 



The terms in curly brackets are the opposite of the Wronskian flHUl) . [11](9.1.16) 

- = -w = — (111) 

TTZo 

and combined with the prefactors yield a —1 which cancels the 1+ in the definition fllOSp of 
the gluon. The final result reads 

z / z ^ 



9e{t>to) = ~TT-\^-j X,{zo)Y^+i{z) (112) 
and has the remarkable property of being factorized in the z- and ;So-dependence. 
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A. 2 £ ^ limit 

In the e —>■ limit of eq. (11121) . both the order and the argument of the Bessel functions grow 
hnearly with t] +oo. By writing the argument in the form 

z = ris, s = 2v^e^*/^ ^ 2^0 , (113) 

using the asymptotic expansions [11] (9.3.6) of the Bessel function J in terms of Airy functions 

Uvs) ^ (t^) ' r ^Ai(r7§C) , iv ^ +oo) (114) 



-(2 = I{s) = / du = I = log VI - s2 , (115) 



Ai(x) ~ _ X 4 exp ( — -x2 ) , (x ^ +oo) , (116) 



and exploiting the asymptotic expansion [11] (10.4.59) of the Airy function 

1 _i / 2 

— p= X 4 exp 

2v^ ^ V 3 

one obtains 

exp [-T] /(2v/^e^*»/2)^ 

•^r^l^o) ~ ^ .1 , (117) 

v/27rr/(l - 4a)4 

where / is the integral defined in eq. 01151) . In the same way, by using the asymptotic expan- 

2 

sion [11] (9.3.6) of the Bessel function Y in terms of the Airy function Bi^rj^Q and the large-i] 
expansion of the latter [11] (10.4.63), one obtains 

« -^^^j^ . (118) 

Before applying the above formulas to the integrated gluon, we can immediately derive the 
£ — i> limit of the unintegrated gluon (cf. eq (lHTl) ) for t > to. In fact 



at > to) = -- (-YYr,{z)Uzo) ^ -f^exp {ry [l{2V^e^'/') - /(2v/^e^W2)] | . 

7] \IJ \Zo J (1 — 4a)2 

(119) 

The difference of the integrals in the exponential yields 



/(2v/Se-/^)-/(2v/Se*'^) = -ir ^^Bl^^ ^ _ y . ,120) 

hence 

n — 1 ^ 



f[t > to) = exp 
V 1 — 4a 



exp[7(a)(t-to)] , (121) 



where 7(a) is the saddle-point value ( l59l) at e = 0. 

As for the integrated gluon density Qs, comparing eq. (11121) with eq. f lll9p we find 

fs{t>to) Z Yr,{z) 



ge{t>to) 2r]Y^+i{z) 



'122) 
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which is nothing but the effective anomalous dimension in eq. (H2l) . We need the asymptotic 
behaviour of 

Y,^^ivs) = Y,{if^ - l)s) = Y,m ^ ^'^^^^ ^^^"^ 



V2^(l -52)4 



where r/ = 77 + 1 and s = s{l — l/rj) ^ s{l — e). From 
f] I{s) ^ (77 + 1) / du = {r] + + Vl - = 77/(5) + log ■ 

Js-es ^ 

and s — > 2^/a, we obtain a finite limit for the ratio (11221) 



(123) 



(124) 



— ~ \^ exp[r] I (s) — 77/(5) 

9e 



a s 



1 - vr^4a 



7(a) (125) 



s=2y/a 



which coincides with the saddle-point 7(a) at e = 0. The explicit expression for the off-shell 
integrated gluon density at e = — namely the gluon in the Qo-scheme — is finally obtained 
dividing eq. fll2ip by 7(a), whence eq. 

In the on-shell case, where the to —00 limit is performed at non- vanishing e, the integrated 
gluon (HTl) at large-?7 behaves like 



9eit) 



a f t 

r exp < — \- n 

7(a)(l-4a)5 12 



!N"(a) exp < i] 



1 



"'da 
2a 



2a 



da 



:i26) 



where I used the Stirling approximation for gamma-functions and eq. (11181) in the asymptotic 
expansion, and an integral representation for the exponent, together with the definition ( l45i) 
for 'N, in the last equality. It is possible to shift the lower limits of integrations to zero, since 
the two logarithmic singularities at a = cancel in the sum of the two integrals. The finite 
additional contribution is provided exactly by the first term "1+" within square brackets. The 
final result is 



Qeit) ^ >J"(a) exp 



1 — -\/l — 4a da 
2 a 



(127) 



Eq. fll27p demonstrates the factorization of the coUinear singularities, and identifies the e-finite 
coefficient factor R = !N(a) and the MS gluon density, according to eq. ( H3l) . 



A. 3 Integrated gluon with running coupling 

In the running coupling case, it is convenient to use the (^-variable introduced in eq. ( l68l) . The 
interesting kinematical region is at t > to 5 which in the regimes Rl and R2 where aB < 1 
corresponds to C < Co < 0: 



Co rC 
jF(x, zo) dx + ^{x, zq) dx 



QeA^ >to) = l- r] 



where the unintegrated density JF is given in eq. (1751) . 

The first integral involves an integral of hypergeometric function of type [11] (15.2.4) 



(128) 



X 



c-2 p I 0-, b 
^ c — 1 



X I dx 



„c-l 



C-1 



2-^1 



a , b 



X 




(129) 
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where the condition of integrabihty at x = is guaranteed by c — 2 = > 0. The second 
integral in eq. fll28p . after the position y = l/x, involves an integral of type [11] (15.2.3) 



y 



a-l 



a + 1 , h 
c 



y 



y F 

a 



a , b 

c 




(130) 



Summing the various contributions yields 

^ 7^(1 -Co) T{2-u)T{l-u) 



9e,b{t > to) 



B T{l-w)T{v -u + 1] 



(131) 



-Co 



w-2 
1 



V - 1 



2-^1 



2-n, 2-v 
3 — w 

2-M, 2-v 
2-w 



Co 2F1 



Co 2-^1 



2F1 



2-u, 2-v 
2-w 



Co ) (-C)'-%i^i 



2 — u , V 

V — u + \ 

2 — u , V — 1 

V — u + 1 

2 — u , V — 1 

V — u + 1 



The hard task is to prove that the (^-independent terms, namely those stemming from the square 
brackets in eq. (11311) . combine themselves in such a way to give a —1 that cancels the 1 at the 
beginning of the r.h.s.. The method is to use relations between contiguous hypergeometric 
functions — differing by one unit in some of their (a, b, c) parameters — and their derivatives. 
By introducing the short-hand notation 



Fi = 2F1 



2-u, 2-v 
2-w 



Cc 



F2 = 2F1 



2 — u , V 
V — u + 1 



1 

Co 



(132) 



and exploiting the relations [11](15.2.6) and [11](15.2.5), one gets 

2-w 



2F1 
2F1 



2-u, 2-v 
3 — w 

2 — u , V — 1 
V — u + 1 



Co 

1 

Co 



{l-u)il-v) 
2-u 



[(i-Co)f;-Fi] 



Coil-u) 



F2 + 



Co-1 
Co^(l-«) 



2 ) 



whence 



(133) 
(134) 

(135) 



In order to find a relation among the F/s and their derivatives, I exploit the Wronskian (J7^ : 



cu 



-Co) 



1—v—w 



1 + V — w 

Co 



F1F2 + F{F2 + -2F1F2 
^0 



(136) 



The combination + FiFl^/Q entering eq. fll35p can thus be expressed in terms of the 

product F1F2. As a result, the terms with F1F2 cancel out and all gamma-functions simplifies. 
Finally, by substituting the explicit expressions (!7T|) of u,v,w, it is straightforward to compute 
the sum of the ^-independent terms in eq. fll3ip and to obtain —1, as I stated previously. The 
remaining (^-dependent term provides the factorized expression (1761) for the integrated gluon 
density. 
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